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ABSTRACT 

A-band galaxy number counts (GNCs) exhibit a slope change at A ^ 17.5 mag not 
present in optical bands. To unveil the nature of this feature, we have derived the 
contribution of different galaxy types to the total A-band GNCs at 0.3 < z < 1.5 by 
redshift bins and compared the results with expectations from several galaxy evolu¬ 
tionary models. We show that the slope change is caused by a sudden swap of the 
galaxy population that numerically dominates the total GNCs (from quiescent E-SO’s 
at A < 17.5mag to blue star-forming discs at fainter magnitudes), and that it is 
associated with a flattening of the contribution of the E-SO’s at 0.6 < z < 1 to the 
total GNCs. We confirm previous studies showing that models in which the bulk of 
massive E-SO’s have evolved passively since z > 2 cannot predict the slope change, 
whereas those imposing a relatively late assembly on them (z < 1.5) can reproduce it. 
The A-band GNCs by redshift bins and morphological types point to a progressively 
definitive build-up of ^50% of this galaxy population at 0.8 < z < 1.5, which can 
be explained only through the major mergers reported by observations. We conclude 
that the slope change in total AT-band GNCs is a vestige of the definitive assembly of 
a substantial fraction of present-day massive E-SO’s at 0.8 < z < 1.5. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: 
formation — galaxies: interactions — galaxies: morphologies 


1 INTRODUCTION 

Galaxy number counts (GNCs, i.e. the number of galax¬ 
ies per apparent magnitude and area unit on the 
sky) were originally conceived to constrain cosmologi¬ 
cal mod els, a purpose for wh i ch they are still being 
used (see Hubble fc TolmarJ 19351; IHartnett fc Hiranolliood 


intermediate-to-high redshifts is still challenging, so most 
studies use p hotometric redshifts to obta i n large complete 


lEnea Romandl201Cli : IShen et al.ll2015l l . However, It was soon 

realized that they are more sensitive to the processes driving 
galaxy evolution than to the assumed cosmo logical model 
ferown fc Tinsle^ll974 iHoversten et al.l[2009f l. 

Although the estimates of luminosity functions (LFs) 
have the advantage over GNCs of showing the intrinsic 
evolution of each galaxy population directly, LF estimates 
are still tricky owing to the assumptions adopted for the 
statistical method, the cosmological model, and the e- 
correction llWillmerlll997ll . as well as the observational se- 
lection effects th at can make any galaxy sample incomplete 
lljohnstonll201ll ~). Completeness in spectroscopic surveys at 
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samples fe.g 

. Perez-Gonzalez et al.ll2005l: iLam et al.l 1201 ll: 

ISobral et al. 

I2OI5IL But small errors in the photometric red- 


shifts translate into large uncertainties in the rest-frame 
fluxes, affecting more noticeably the LF estimates when the 
object statistics are low. Consequently, the LFs of bright (i.e. 
massive) galaxies tend to be poorly estima t ed and their evo¬ 
lution is less constrained llGoto et al.ll201ol : iKhostovan et all 
I2OI5II . Additionally, there are very few studies deriving k- 
corrections from real data, and they are usually limited to 
low-to-intermediate redshifts, certain bands, and/or some 
specific galaxy t ypes, so LFs are affected by these uncer¬ 
taint i es too (e.g. Rocl2eetjjj [2009l : iHuertas-Companv et al.l 


I2OO9I : IChilingarian et alT 2010l L These uncertainties in the 
fc-corrections also increase the difficulty of modelling. In 
fact, current models can predict the galaxy luminosity func- 
tions at z ^ 0, b ut not their precise evolution with redshift 
iIGuo et al.ll201lh . 

On the contrary, GNCs are free from these uncer¬ 
tainties and have higher statistical significance, so they 
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are quite robust observational diagnostics for discriminat¬ 
ing between different evolutionary models.Nevertheless, de¬ 
spite being apparently simple, GNCs summarize the ef¬ 
fects of different galaxy evolutionary processes taking place 
at different redshifts, so they still represent a challenge 
for current galaxy evolution models, which have difficulty 
in rep roducing total GNCs from NUV to NIR simultane - 
ously dKitzbichler fc Whitell2006l . l2007l : [Grazian et aLl2009ll . 
This difficulty arises because the slope of GNCs in op¬ 
tical bands is almost featureless and nearly Euclidean 
for a wide magnitude range ( d\o^(N)/dm ~ 0.6, see 
[Peebles! I I 971 I : lYasuda et al.l[200lh . but it exhibits a notice- 
able slope change i n NIR GNCs at K 17.5 mag (see 
Gardne£_et_^ ]ll99^jBershudj_et_^lJ[w9^lMcCrackene£^lJ 


2000|^ Cristbbffi-H ornillos et al.ll200S . l2009| -1l mai et al. I I 2 OO 7 I: 

Barro et al. 20091 7 This slope change is also reported 


in the J and H bands^ alth ough it is less noticeable 
lICristobal-Hornillos et ^120091 1. 

Many studies have tried to derive the processes re¬ 
sponsible for this feature through modelling, but with lit¬ 
tle success. In general, whereas UV and blue optical GNCs 
require strong luminosity evolution of g alaxy populations 
and/or high merger rates to be fitted (lHogg_et_alJ_[l993 
Gardner_e£alJ_ 20001: [Fried et ahl I 2001 I : iNagashima et cil] 


2 OO 2 I : iioversten et al.l 2009fl . r^ optical and NIR GNCs 


basically follow a pure luminosity evolution (PLE) trend 


slope chanEe. see Pozzetti et al.l 19961: He & Zhand 

1998 

Totani & Yoshiil I 2 OOOI: Martini 

I 2 OOII: iHuane: et al.l 

2001 

Martinil 2001: iTotani et al. 2001 

ISomerville et al.ll200' 

). No 


consistent evolutionary scenario has so far been proposed to 
reproduce GNCs simultaneously by redshift bin and mor¬ 
phological type in all bands. 

iBershadv et al.l lll998!) compared deep GNCs in the J 
and K bands with the expectations of different galaxy evo¬ 
lutionary models. These authors had already attributed the 
slope change in the NIR bands to the evolution experienced 
by the ancestors of present-day bright quiescent galaxies, be¬ 
cause they had found a deficit of galaxies with J — K colours 
typical of bright passively-evolving objects at 1 < 2 < 3 with 
respect to non-evolutionary models that indicated that these 
gala xies were much bluer a t these epochs than at 2 < 1 (see 
also lGardner|[l999ll . It was ICristobal-Hornillos et al.l ll2003l . 
CH03 hereafter) who first suggested that the slope change 
could be reproduced by delaying the formation of massive 
early-type galaxies (E-SO’s with masses M* > 10^^M©) 
down to 2 < 2. However, their model failed to reproduce 
GNCs in blue optical bands unless a transient population of 
star-forming galaxies was included ad hoc at all redshifts. 
In any case, the assembly epoch proposed for massive early- 
type galaxies was strikingly late compared to paradigms 
(mainly monolithic collapse at 2 > 3) that were prevalent in 
those years. 

iKitzbichler fc Whitel 1 I 2 OO 6 !) showed that PLE models 
for early-type galaxies overpredict the population of massive 
E-SO’s at high redshift. They supported the conclusions by 
previous authors that, in order to reproduce GNCs data, ei¬ 
ther the majority of massive E-SO’s should be assembled at 
late epochs, or the ancestors of these galaxies should have 
suffered strong dust extinction phases at 2 > 1. A coeval 
study proved that both mechanisms were indeed required 
to reproduce total GNCs simultaneously in the U, B, and 


K bands llEliche-Moral et al.ll2006H EM06 henceforth). The 
appearance of the bulk of massive E-SO’s on the cosmic sce¬ 
nario at a quite late redshift (2 = 1.5) was necessary to 
reproduce the slope change in the NIR GNCs, but a moder¬ 
ate dust extinction in the early evolutionary stages of these 
galaxies was necessary to remove ‘bump’ signatures in the 
GNCs of blue bands caused by the young stellar populations 
contained in these galaxies at 2 ~ 1.5. The EM06 model did 
not identify the mechanism responsible for the E-SO assem¬ 
bly, but these authors suggested that the high extinction 
required in the early phases of this galaxy population may 
indicate gas-rich major mergers. 


Later studies based on previous models estimated that 
the present-day number density of massive early-type galax¬ 
ies should decrease by a factor of -^2-3 at 2 ^ 1 to repro¬ 
duce multiwavelength GNC data (ICristobal-Hornillos et aP 
I 2 OO 9 II . This association of the slope change in NIR GNCs and 
the late bui l d-up o f massive red objects were supported by 
I Barro et al.l (l2009t) . who showed via GNCs by redshift bins, 
that the slope change in NIR GNCs was caused by a promi¬ 
nent decrease in the number density of L ~ L* objects at 
0.8 < 2 < 1. However, neither the galaxy population nor the 
mechanism responsible for the slope change in the A'-band 
GNCs have been observationally identified or conhrmed yet. 


In order to test observationally whether the slope 
change in NIR GNCs is really a trace of the recent build-up 
of massive E-SO’s or not, we have analysed the contribution 
to the total A'-band GNCs of the various galaxy morpholog¬ 
ical types (in particular, the E-SO’s) at 0.3 < 2 < 1.5. The 
results obtained have been compared with the predictions 
of several galaxy evolution models in order to shed some 
light o n the mechanisms resp onsible for this feature in NIR 
GNCs. iGrazian et al.l (I 2009 I 1 proposed that the discrepan¬ 
cies between models and GNC data could be evaluated in 
more de tail by dividing the counts into redshift bins, as was 
done bv iBarro et al.l (l2009l L GNCs in redshift bins can dis¬ 
entangle global evolutionary effects that may depend on the 
redshift, and that are mixed in the total GNC, thereby at the 
same time avoiding the large uncertainties associated with 
LF estimates (see above). We have therefore also derived the 
GNCs of E-SO’s by wide redshift bins in our sample in order 
to compare with the expectations from the models. 


The paper is organized as follows. In )J2]we describe the 
data and the procedure used to derive GNCs by redshift 
bin and by morphology in the K band. In jjS] we provide a 
brief outline of the GNC models with which we have com¬ 
pared our observational results. The results are described 
and discussed in El Section |4T] discusses the results de¬ 
rived directly from the observational GNCs obtained here 
concerning the nature of the slope change in NIR GNCs. 
In il4.21 we show the clues to the evolutionary mechanisms 
of massive E-SO’s that arise from the comparison of GNC 
data with the predictions of the models. The conclusions 
are hnally provided in !j5] We use a concordant cosmology 
throughout the paper (Hm = 0.3, H a = 0.7, Hq = 70 km 
s“^ Mpc“^, see ISoergel et al.|l2007^ . All magnitudes have 
been converted into the Vega system. 
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Figure 1. Colour-magnitude diagrams for the P13 observational sample with K < 19.5 mag in the three wide redshift bins (0.3 < z < 0.7, 
0.7 < 2 : < 1.1, and 1.1 < 2 < 1.5) considered in that study. Red and blue galaxies, according to the criteria used by P13, are plotted with 
red crosses and blue diamonds, respectively. The median errors for the observational data in each axis are plotted in the bottom-left 
corner of each panel. Top panels: rest-frame (U — B) colour vs. rest-frame absolute magnitude in the B band, Mg^rest—frame? for th® 
galaxies in the P13 sample, using t he values available i n the Rainbow database. The red dotted line corresponds to the fit to the red 
sequence E-SO galaxies obtained by Ivan Dokkum et al.l ll200d') , and the green solid line is the same line shifted downwards by 0.25 mag 
to p ass through the vall ey between the red and blue galaxies (usually used as a cut to distinguish between red and blue populations; 
e.g. IWillmer et al.|[2006l^ . Bottom panels: available apparent colours nearest to the rest-frame (U — B) colour at each redshift interval 
vs. Mg rest—frame for the Same galaxies. Horizontal red solid lines: colour cuts defined by P13 to distinguish between red and blue 
galaxies on the basis of apparent colours. Lines limited by stars: theoretical trends followed by different galaxy types within the 
redshift interval of each frame, modelled assuming standard SFHs starting at 2 / = 3 and char acteristic physical propert ies according to 
observations. These trends have been obtained with the stellar population synthesis models bv iBruzual Sc Char Ion ll2003l , see the text for 
details). A stellar mass of log(M*/MQ) = 11 has been assumed for all models. Red: E—SO galaxy. Yellow: Sa-Sab. Green: Sb-Sbc. Blue: 
Sc-Scd. Purple: Sd-Irr. Pink: Dust-reddened starburst galaxy with Ay = 3.0 mag. [A colour version is available in the online version.] 


2 DATA AND GNC ESTIMATION 

This study is based on the data, s election, and classification 
procedures by iPrieto et alJ (|2013l . P13 hereafter). Detailed 
information is provided in the original paper, so here we 
provide only a brief description. 

We have used broad band photometric data {U, B, V, 
J, J, and K) from the Rainbow Extragalactic Databasqj 
(|Barro et al.l 120113) over an area of ~155arcmin^ in the 

^ Rainbow Extragalactic Database: 

https: //rainbowx.fis.ucm.es/Rainbow_Database/Home.html 


Groth Strip (a = 14'"16'”38.8° and 5 = 5 2° 16'52", see 
Groth et al.l Il994l : iRatnatunga et al] 1 19991 : ISimard et al.l 
20021) . The database contains photometric data from the 
UV to the FIR over this field, as well as diagnostics derived 
from the fit of the spectral energy distributions of nearly 
80 000 galaxies selected in IRAG bands. We have nsed pho¬ 
tometric redshifts available in this database, with an accu¬ 
racy of (A2/(1 -|-z)) ~ 0.03 (see fig. 2 in P13). The precision 
of the photometric redshifts (measured as the 68% confi¬ 
dence interval in the probability distribution function of the 
redshift estimates around the most probable photometric 
redshift computed with the Rainbow code) differ according 
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to reds hift, but are typic ally below Azphot 0.08 (see ta¬ 
ble 2 in iBarro et al.llioilbll . The depths reached are typically 
B ~ 26 and K 21 mag. We started with a sample of ~850 
sources with K < 19.5 mag at 90% completeness in terms of 
detection. This sample is complete in detection at the 100% 
level for K < 18.5 mag, thus ensuring total completeness of 
the GNCs down to this limiting magnitude. Consequently, 
the data adequately cover the range of magnitudes in which 
the slope change in the Tt'-band GNCs appears without re¬ 
quiring efficiency corrections. We preferred to avoid such 
corrections because they often fail to recover the true GNC 
distribution owing to data errors and unrealistic assump¬ 


tions adopted in the models develo ped to estimate them (see 
comments concerning this in, e. g., ICristobal-Hornillos et al.l 
l2008l : lEliche-Moral et al.ll2006bh . 

In order to understand the differences between the 
GNCs in NIR and UV-blue bands, we first classified the 
galaxies according to their colours into those belonging to 
the red sequence or the blue cloud. This bimodality in the 
basic properties of galaxies can be traced up to z ^ 2 
llKodama et al.ll2007l: I Williams et al.ll2009l : iBrammer et al.l 
I 2 OIII : Whitaker et al.ll2011^ . The rest-frame U — B colour is 
usually used for this purpose. This colour traces the 4000 A 
break, a spectral feature characteristic of evolved stellar pop¬ 
ulations, so it is a good discriminator between really quies¬ 
cent (red) and star-forming (blue) galaxies. However, we had 
to perform this classification in terms of apparent, instead 
of rest-frame, colours because of the large uncertainties in 
these latter that result from the propagation of the errors in 
the photometric redshifts and magnitudes. 

We illustrate this problem in Figure [T] In the top pan¬ 
els, we plot the distribution of our galaxies in the rest- 
frame U — B vs. Mb diagram for three wide redshift bins 
(0.3 < z < 0.7, 0.7 < z < 1.1, and 1.1 < z < 1.5). We 
have overplotted the linear fit performed to descr i be th e lo¬ 
cation of the red sequen ce by van Dokku m et al.l (I 2 OOC 1 I I. as 
well as the limit used bv lWil'lmer et al.l ( 2006ll to isolate the 
galaxies on the red sequence from those lying in the blue 
cloud. We have represented the median error bars in both 
axes in each panel. Attending to the typical errors in our 
rest-frame U — B estimates, it is obvious that we cannot 
robustly assign a colour class to the galaxies on the basis 
of a colou r-magnitude cu t in this diagram as in the one 
defined by IWillmer et al.l . This prevented us from trusting 
these rest-frame colours for making the red-blue distinction 
in our sample. 

We used instead the observed apparent colour from our 
dataset that traces the rest-frame U — B more closely at the 
centre of the three wide redshift bins under consideration 
(U — V, B — I, and V — K for the low-, middle-, and high- 
redshift intervals, respectively). We show the distribution of 
these colours for our galaxies as a function of the rest-frame 
Mb in the bottom panels of Fig.[T] for the three redshift 
bins. The apparent colour distributions also exhibit the well- 
known bimodal distribution of galaxies into red and blue 
populations (as do the rest-frame U — B diagrams), with 
the advantage of having much lower errors in the colours 
(compare the median error bars in these diagrams with those 
in the corresponding rest-frame diagrams at the top panels). 
Therefore, the classification into red and blue galaxies is 
much more reliable for our data if it is defined according to 
these apparent colours. 


We thus defined the colour cuts assuming Gaussian- 
shaped distributions for the red and blue galaxy popula¬ 
tions in the histogram of the apparent colour corresponding 
to each redshift bin. The cut was defined as the point where 
the fitted Gaussian distributions crossed. The blue and red 
distributions were clearly defined in the first two redshift 
bins, but not in the last one. We defined the colour cut 
in the highest redshift interval by ensuring that it selected 
similar galaxy types to those defined in the lower redshift 
bins. For this purpose, we used the evolutionary tracks fol¬ 
lowed in colour-redshift diagrams by characteristic galaxy 
types, as derived from the s tellar population synthesis mod¬ 
els by iBruzual &: Gharlotl (l2003li . The resulting cuts were 
{U — B) = 1.4 mag for 0.3 < z < 0.7, {B — I) = 2.7 mag for 
0.7 < z < 1.1, and {V — K) = 4.5 mag (see Fig. 3 in P13). 

In Fig.m we have also plotted the locations expected 
for some of these models in the colour-magnitude diagrams 
for each redshift interval. We overplot the colours and mag¬ 
nitudes expected for a typical E-SO, Sa-Sab, Sb-Sbc, Sc- 
Scd, Sd-Irr, and dust-reddened starburst galaxies, as de¬ 
fined by the parameterizations provided for the star forma¬ 
tion histories (SFHs) and the characteristic physical proper¬ 
ties (metallicity, dust extinction, e-foldi ng timescale for the 
star for mation) of each type defined in lEliche-Moral et all 
(l2010bl . EMIO henceforth), based on observational studies 
(see S 43.1l and l3.4l for more information). The star formation 
is set to start at z/ = 3 for all types, and a stellar mass of 
10^^ M© is assumed for them all. The top panels of the figure 
shows that E-SO’s and dust-reddened star-forming galaxies 
are expected to lie in the region of ‘red galaxies’ according 
to the colour cuts defined in the rest-frame U — B diagrams 
for all redshifts, whereas early-type spirals would lie in the 
green valley and late-ty pe spirals in the blue cloud. Accord¬ 
ing to the cut defined bv lWillmer et all ll2006li . E-SO’s, dust- 
reddened starbursts, and Sa-Sab’s would be identified as 
‘red galaxies’, and spirals of types later than Sb lie would be 
‘blue galaxies’. In the bottom panels, the equivalent colour 
cuts defined in P13 in the apparent colours basically distin¬ 
guish between the same galaxy models, thus ensuring that 
these colour cuts select homogeneous samples of galaxy pop¬ 
ulations in the whole redshift interval in a similar way to 
those defined in terms of rest-frame colours (see Figs. 3 and 
4 and §3.1 in P13 for more details). 

Additionally, colour-colour diagrams were used to re¬ 
move stars from the sample. We have finally got a red galaxy 
sample with 302 systems and a blue one of 536 objects with 
K < 19.5 mag. The blue galaxy sample basically consists 
of star-forming spirals. The P13 study was centered on the 
build-up of the red sequence, so no morphological classifica¬ 
tion was performed for blue galaxies. In the present study, 
we show that the slope change in the K band is not related 
to this galaxy population, in the sense that the GNGs of blue 
galaxies do not exhibit any change of trend at 16 < A" < 19.5 
(see 1)4.111 . Therefore, we will consider that blue galaxies 
make up an independent galaxy type throughout the paper 
and we will not distinguish different morphological types 
among them. 

In the red sample, however, we have used the visual 
morphological classification developed by P13. The classifi¬ 
cation criteria were centered on the global morphology of 
the galaxy (if it was a spheroid, or a disc-dominated object) 
and the distortion level of the whole galaxy body (whether 
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the galaxy was ‘regular’, in the sense of poorly distorted 
morphology, or ‘irregular’, meaning that it exhibited a sig¬ 
nificant distortion). The definitive visual classification ac¬ 
counted for the median results of three independent classi¬ 
fiers. We used surface brightness isophotes and surface maps 
in the I band (besides the images) to identify structural dis¬ 
tortions and morphological features. This band was selected 
for the visual classification because it is the reddest optical 
band with the best spatial resolution available in the data 
set. 

Spheroid-dominated galaxies were defined as those ex¬ 
hibiting central light concentrations resembling those of 
nearby E-SO’s (i.e. the central bulge controlled the surface 
brightness profile up to at least one third of the galaxy size). 
We identified as ‘irregular or distorted’ galaxies only those 
with morphological features characteristic of major mergers, 
such as train-wreck morphologies, equal-length tidal tails, 
or double nuclei objects with high and similar b rightness, 
following the methodology by Ijoeee et al.l ll2009ll . We thus 
avoided the inclusion of galaxies at an early stage of a ma¬ 
jor encounter, undergoing minor mergers, or experiencing 
inhomogeneous starbursts along their discs into this ‘irreg¬ 
ular’ class. However, some contamination by these objects 
cannot be excluded. Consequently, attending to the global 
morphology and structural distortion level, P13 defined six 
major exclusive classes of red galaxies: 

• Compact galaxies.- Galaxies exhibiting compact mor¬ 
phologies, according to the seeing of the images. The number 
of these objects in the sample was negligible. 

• Regular Spheroids.- Galaxies dominated by a central 
spheroid, with regular isophotes. 

• Regular Discs.- Galaxies dominated by a disc compo¬ 
nent with regular isophotes. 

• Irregular Spheroids.- Galaxies with irregular isophotes, 
dominated by an spheroidal component. 

• Irregular Discs.- Galaxies with irregular isophotes, 
dominated by a disc. 

• Non-Classified.- Galaxies that could not be classified 
owing to their faintness or noise. The number of these ob¬ 
jects was negligible too. 

The visual classification was validated by comparing the 
results with the automatic classifications resulting from two 
quantitative parameters known to be proxies of the concen¬ 
tration and distortion level. The tests developed to estimate 
the robustness of the visual classification are detailed in §4.2 
of P13. 

We have estimated the GNGs for red regular spheroid- 
dominated (basically, E-SO’s), and red regular disc- 
dominated (hereafter, ‘red disc’) galaxies, as well as for red 
galaxies with highly distorted morphologies and blue ones. 
The GNGs obtained at 16 < A" < 18.5 mag are listed in 
Table[TJ both in total and for each galaxy class. The er¬ 
rors listed in the Table account quadratically for the statis- 
tical counting uncertainties at the 84.13% confidence level 
llGehreIslll986 1 and the classification errors. The uncertain¬ 
ties due to classification are estimated as the square root 
of the quadratic mean of the statistical errors of the three 
independent classifications. 

In Table[5] we list the /f-band GNGs that we have 
obtained for red regular spheroidal-dominated systems (E- 
SO’s) in different redshift bins and for the whole redshift 


range from 2 = 0.3 to 2 = 1.5. The errors listed in the ta¬ 
ble correspond to the quadratic summation of the squared 
uncertainties associated with statistical counting, classifica¬ 
tion error, and uncertainties associated with cosmic vari¬ 
ance, which have been estimated as described below. We 
have not accounted for the uncertainties due to redshift er¬ 
rors because they were negligible for the redshift bins under 
consideration (see P13 for more details). 

We have derived cosmic variance estimates associate d 
with our sample using the model by iMoster et al.l ll201lfl , 
which predicts the cosmic variance associated with a galaxy 
population of a given stellar mass considering a ACDM cos¬ 
mology and the galaxy bias. These authors have kindly made 
a cosmic variance calcula tor publicly available (QUICKCV, see 
iNewman fc Mosterj|2014l ~) , which provides the root of the cos¬ 
mic variance of galaxies (ugg) expected for different stellar 
mass bins in different redshift intervals for a survey of a 
given geometry. 

The objects dominating the GNGs of quiescent galax¬ 
ies in each redshift bin are those with L ~ L*, because the 
LF of these objects always has a decreasing slope at faint 
magnitudes at all redshifts up to 2 ~ 2 (in fact, a > —0.6 
at 2 < 2 for this galax y population, see lArnouts et aDl2007l : 
ICirasuolo et al.ll200l1l . The stellar mass typical of quiescent 
galaxies with L ^ L* rises from ~ 2 x IO^'^Mq to ~ 3- 
5 X 10^^M© from 2 0 to 2 ~ 1.5, considering the typical 

brightening experienced by this galaxy population and the 
evolution of the mass-to-light ratio in the K band derived for 
these objects bv lArnouts et af] ll2007lf . This means that the 
objects that control the GNGs of E-SO’s at all magnitudes 
(L > L*) have stellar masses M, > up to 2 ~ 1.5. 

For galaxies with stellar masses log( M«/Mf?)) > 11 (hence¬ 
forth, simply ‘massive’), the model bv iMoster et al.l predicts 
that the fractional error is up to >^28-30% for the 

redshift bins and the field considered here. Therefore, our 
GNGs of E-SO’s by redshift bin are affected by cosmic vari¬ 
ance uncertainties of ~30% at most. We have included this 
maximum uncertainty in the errors of the GNGs of E-SO’s 
by redshift bin in Table[21 as noted above. 

Tabled also shows results on the GNGs of E-SO’s by 
redshift bin, but for wider redshift intervals (0.4 < 2 < 0.8, 
0.6 < 2 < 1.0, and 0.8 < z < 1.2) than those used in Ta- 
ble[21 We have derived them in order to reduce the errors 
associated with the distributions of GNGs of E-SO’s in nar¬ 
row redshift bins, by raising the statistics and reducing the 
cosmic variance uncertainties. This allows us to derive more 
robust conclusions from the comparison with the models, 
as we show in 114.2.31 The errors in Table[3] also account for 
statistical counting and classification errors, as well as for 
cosmic variance uncertainties, which are reduced to <15% 
for these redshift bins. Note that these wide redshift inter¬ 
vals overlap. 


3 GNC MODELS 

We have compared the results obtained on GNGs in the 
K band in Sect .[2] with the predictions of different evolu¬ 
tionary models, in order to shed some light on the evo¬ 
lutionary mechanisms responsible for the slope change at 
K ~ 17.5 mag. The models used for comparison were the 
following: 
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Table 1. iC-band GNCs, both total and by galaxy type as defined in P13 


Red normal bulge-dominated* Red highly distorted Red normal disc-dominated Blue galaxies Total 


K* 

(1) 

log 

(2) 

(3) 

aJ 

(4) 

log N 

(5) 

Ai 

(6) 

Au 

(7) 

log N 
(8) 

Ai 

(9) 

Au 

(10) 

log JV 
(11) 

Ai 

(12) 

Au 

(13) 

log JV 
(14) 

Ai 

(15) 

Au 

(16) 

16.25 

2.79 

0.09 

0.09 

1.80 

0.30 

0.44 

2.11 

0.21 

0.26 

2.58 

0.18 

0.12 

3.08 

0.10 

0.08 

16.75 

2.82 

0.09 

0.09 

2.11 

0.21 

0.26 

2.28 

0.17 

0.20 

2.76 

0.15 

0.11 

3.19 

0.09 

0.07 

17.25 

3.15 

0.06 

0.06 

2.65 

0.11 

0.12 

2.58 

0.12 

0.13 

3.19 

0.09 

0.07 

3.58 

0.06 

0.05 

17.75 

3.21 

0.06 

0.06 

2.98 

0.07 

0.07 

2.50 

0.13 

0.14 

3.42 

0.07 

0.06 

3.74 

0.05 

0.04 

18.25 

3.11 

0.06 

0.06 

3.25 

0.05 

0.05 

2.65 

0.11 

0.12 

3.66 

0.05 

0.05 

3.91 

0.04 

0.04 

18.75'’ 

3.09 

0.06 

0.08 

3.44 

0.04 

0.04 

- 

- 

- 

3.83 

0.04 

0.04 

4.03 

0.03 

0.03 


* For the equivalence of the classes defined in P13 with traditional Hubble types, see m 

* X-band apparent magnitudes in the Vega system. 

t The GNCs are provided in units of mag“^ deg“^. Their logarithmic values are listed. 

t The A; and Au errors are defined as A; = log(iV) — log(7V — cti) and Au = log(7V + CTu) — log(Al), where (Tj and CTu correspond to the 
lower and upper errors due to statistical counting and classification uncertainties in each case. 

^ The data in this last magnitude bin are not complete in detection. 


(i) a mo del assuming: PLE for all g alaxy types, including 
the E-SO’ s (lEliche-Moral et ^ 12010a); 

(ii) the iKitzbichler fc White (2007) model (KW07 here¬ 
after), based on the ACDM scenario of the Millenium simu¬ 
lation; 

(iii) the EM06 model, assuming a recent buildup of E- 
SO’s at 2 = 1.5 and moderate dust extinction at the initial 
evolutionary stages of these galaxies; and 

(iv) the EMIO model, proposing a hierarchical scenario 
for the formation of massive E-SO’s based on observational 
major merger rates up to z 1.5. 


The two first ones are already kn own to not reproduce 
the s l ope ch ange of the Tf-band GNCs (IKitzbichler fc Whitel 
l2006l . [ioOTh . but they are used here as a reference. We have 
compared the expectations of these four basic models with 
our observational GNCs in m Below, we provide a brief de¬ 
scription of the basic assumptions and main characteristics 
of each model. 


3.1 PLE model for E—SO’s since 2 = 3 


The PLE model assumes that the stellar populations of all 
galaxy types have just evolved according to an SFH char- 
acteri stic for each type, but no number evolution is consid¬ 
ered dEliche-M oral et al.l l201()i] l. We have used the NCMOD 
code by iGardneil ( IQQsf l to predict the GNCs in different 
bands. The code starts from the local LFs of different galaxy 
types and evolves them backwards in time by considering 
three sources of evolution: the number evolution according 
to Gardner’s implementation of merging (it is switched off in 
this case of PLE), the typical luminosity evolution of each 
galaxy type due to the SFH assumed for it (including k- 
and e-correct ions), and the evolution of the volume element 
derived from the assumed cosmology. The ‘counter-evolved’ 
LF of a given morphological type at each redshift is used 
to estimate the contribution of these galaxies to the total 
GNGs at a given magnitude and photometric band. 

The galaxy types considered at 2 ~ 0 in the PLE model 
have been the same as those used in EMIO (E-SO, Sa-Sab, 
Sb-Sbc, Sc-Scd, and Sd-Irr). The Schechter functions of the 
LFs at 2 = 0 for each type are described in table 1 of EMIO. 
The SFH of each galaxy type is fixed through the star for¬ 
mation rate (SFR) assigned to each type and the redshift at 


which the star formation activity in each galaxy class starts 
(zf). In general, exponential functions have been used to 
represent the SFRs of spirals, instantaneous bursts for E- 
SO’s, and constant SFRs for Sd-Irr and Im types. Different 
exponentially decaying timescales have been considered for 
spirals, according to observations (r = 4, 7, and 9 Gyr for 
Sa-S ab, Sb-Sbc, and Sc-Scd respectively, see lGallazzi et al.l 
I 2 OO 5 II . The characteris tic metallicity of each type has also 
been fixed according to lGallazzi et al.l (l2005li . The complete 
set of parameters for describing the SFH and metallicity of 
each type can be found in table 2 of EMIO. The dust extinc¬ 
tion for each type was set to typical observational values, 
as in the EMIO model (see EH for more details). We have 
assumed 2 / = 3 for all types. 


3.2 Hierarchical assembly of E—SO’s based on the 
Millenium Simulation (KW07 model) 

In order to test a hierarchi cal scenario of galaxy formation, 
IKitzbichler &: Whitel (l2007l . KW07 hereafter) compared the 
GNCs predicted in the ACDM cosmological framework of 
the Millenium simulatior0 with data in various bands. They 
derived the GNCs by constructing deep light-cone surveys 
from the attaching of semianalytic models of galaxy for¬ 
mation to the merging trees of the dark matter derived in 
the simulation (ISoringel et al.ll2005h . The models accounted 
for the baryonic processes in the haloes and sub-haloes, in¬ 
cluding the ‘radio mode’ feedback from the central galaxies 
in groups and clusters, and were very successful in repro- 
ducing several global properties of the galax y population 
dCroton et al.ll2C)(5^ : IPe Lucia &: BlaizHl2007^ . 

KW07 already showed that this model reproduced the 
data of GNCs in optical bands, but predicted an overabun¬ 
dance of massive galaxies at 2 > 1 that were too red, 
despite the hierarc hical nature of the scenario (see also 
iGrazian et al.l I 2 OO 9 I L We have compared our results with 
the total GNCs in several bands derived by these authors. 


^ Millenium Simulation Project webpage available at: 
http: / / WWW. mpa-garching. mpg.de / galform / virgo / millennium / 
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3.3 Ad hoc buildup of E—SO’s at z = 1.5 (EM06 
model) 

The EM06 model considered that the whole population of 
early-type galaxies appears on the cosmic scene at the same 
epoch {z 1.5), but the model does not assume any mecha¬ 
nism for it. This framework is an obvious over-simplification 
of reality because the history of the specific SFR of present- 
day E-SO’s points to a progressive build-up of their stellar 
populations, with sta r-formation peaks at some epochs de - 
pending on the mass (iThomas et al.l[2010l : IChoi et al.ll2014ll . 
However, this model mimics a late structural assembly of 
E-SO’s by approximating the redshift at which this build-up 
takes place by the redshift in which the bulk of their stars 
are born (zf = 1.5 in this model). 

The model is based on the original implementation of 
the NCMOD code (see for a brief description). It as¬ 

sumes four galaxy types (E-SO, SO/a-Sb, Sbc-Sd, and Im), 
taken from the l ocal L Fs derived using SDSS data from 
iNakamura et al.l (120031 '). Note that these are coarser than 
those used in the PLE and EMIO models (their Schechter 
parameterizations are fisted in table 11 of EM06). The red¬ 
shift formation of all spiral types has been set to Zf = 3, 
whereas Zf = 1.5 for E-SO’s. Similar SFHs and metallicities 
to those commented on in the PLE model were adopted in 
this model, including instantaneous starbursts for the E-SO 
type (consult table 12 in EM06). 

The merging procedure implemented in the original 
NCMOD code was assumed but only affected the faint end 
of the GNCs (and that very negligibly), so we can ne¬ 
glect their effects here (see comments in EM06). The effects 
of dust extinction were already considered in the original 
NCMOD code. An average value of the B-band optical depth 
of Tdust.star = 0.6 was assumed for all morphological types 
(also for E-SO’s). This value implies the consideration of 
strong dust extinction in the early phases of E-SO’s, which 
was essential to avoid the appearance of a bump at faint 
magnitudes in the GNCs of blue optical bands caused by the 
blueish colours of these galaxies during the ~lGyr period 
after their build-up at 2 ; = 1.5 (see Q. As the EM06 model 
is not limited in magnitude range, it provides estimates on 
GNCs for all apparent magnitudes. 


3.4 Progressive hierarchical assembly of E—SO’s 
(EMIO model) 

In order to test the feasibility of the hierarchical origin of 
massive E-SO’s (M* > 10^^ Mq), EMIO studied how the 
population of massive early-type galaxies would evolve back¬ 
wards in time under the assumption that they have derived 
from the major mergers that are obser vationally reported 
at eac h redshift up to 2 ~ 1.5 (see also lEIiche-Moral et al.l 
l2010d ). Therefore, the EMIO model simulates a more re¬ 
alistic scenario for the build-up of massive E-SO’s than 
the EM06 model, and distinguishes between the redshift of 
structural assembly of an E-SO (2asaembiy) and the redshift 
at which the bulk of stars in the final galaxy formed in its 
progenitors ( 2 /). In this model, the build-up of the local 
massive E-SO population is progressive and can span a wide 
range of redshift values 2assembiy. 

The basic assumption of the model is that each major 
merger registered by observations has given rise to an E-SO 


galaxy, which is a reasonable hypothesis accounting for the 
typical gas contents of galaxies at 2 < 1.5 (see references 
in EMIO). The model accounts for the relative contribution 
of dry, mixed, and wet mergers, the existence of transient 
stages of major mergers with strong dust extinction and 
irregular morphologies (in which both galaxies have already 
merged into one body), and the typical timescales of these 
phases. 

The model predictions are limited to Mb < —20 mag, 
which is the limiting magnitude of the observational sample 
used to derive galaxy merger rates. This limiting magni¬ 
tude implies that the model can only trace back in time the 
evolution of E-SO’s brighter than this magnitude at all red- 
shifts because this is the only galaxy type assumed to be 
affected by the ‘inverse merging’ of the model (i.e. it is the 
only one assumed to derive from major merging). This lim¬ 
its the model predictions to the E-SO population that ends 
up with > IO^^Mq at 2 = 0 (see EMIO for more de¬ 
tails). Therefore, the E-SO’s sampled by the model at each 
redshift present lower masses as the model simulates the 
evolution backwards in time, not just because of the inverse 
SFH assumed for them (which is passive), but also because 
they include the gas-poor progenitors of the decomposed 
E-SO’s, which obviously have lower masses than the rem¬ 
nant E-SO. The mass ratio of each encounter is randomly 
distributed between 1:1 and 3:1, meaning that the masses 
of the progenitors range from ~30% to ~60% of the stel¬ 
lar mass of the E-SO resulting from the merger. The model 
distributes the estimated total number of major mergers oc¬ 
curring at each redshift according to the trends with the 
mass derived from observation s and cosmological si mula- 
tions dPe Lucia fc Blaiz^l2007l : iFontanot et al.ll2009ll . The 
remaining Hubble types under consideration experience only 
luminosity evolution, except for the late-type discs, which 
also increase its population at the expense of the gas-rich 
progenitors of the E-SO’s that are being decomposed at each 
redshift. 

The model was computed with an improved version of 
the NCMOD code. The local galaxy types considered were thin¬ 
ner than those assumed in the EM06 model (E-SO, Sa-Sab, 
Sb-Sbc, Sc-Scd, and Sd-Irr). The parameterizations of the 
LFs, SFHs, and metallicities of each type are the same as 
those described in the PLE model (asm- All types have 
the redshift of start of their SFHs at 2 / = 3 (including the 
E-SO’s). 

In this model,which is strictly based on observations 
and on robust computational results, a more realistic treat¬ 
ment of the effects of major mergers on the LF of E-SO’s 
was implemented. The model considered the different phases 
and time-scales in a major merger. Because the structural 
distortion and dust extinction of a merger remnant depend 
strongly on the characteristics of the progenitor galaxies, 
the model considered different properties of the E-SO rem¬ 
nants during their transient relaxation phase, depending on 
whether they derive from a gas-rich or gas-poor encounters, 
according to the observed f raction of each merger type tak¬ 
ing place at each redshift (iLin et al.ll2008ll . The model as¬ 
sumed that the galaxies undergoing transient intermediate- 
to-advanced stages of major mergers had globally distorted 
morphologies and should be very red because of the strong 
dust reddening associated with the merger-induced star- 
bursts in gas-rich events, or owing to the lack of star for¬ 
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mation in gas-poor mergers. This is why the galaxies pre¬ 
dicted to be undertaking an advanced stage of a gas-rich ma¬ 
jor merger were called dust-reddened star-forming galaxies 
(DSFs) in EMIO. However, because we have dust-reddened 
starburst galaxies in our sample that are not necessarily re¬ 
lated to major mergers, we will refer more specifically in the 
model to these galaxies as advanced stages of major mergers. 

The dust extinction procedure in the original NCMOD 
code was improved by allowing the various galaxy types to 
acquire a different Tdust.star according to observational values 
(see table2 in EMIO). The E-SO’s were assumed to have lit¬ 
tle dust extinction, in better accordance with observations 
("Tdust.star = 0.1 in the B band). For the intermediate-to- 
advanced stages of gas-rich major mergers, the model as¬ 
sumed a high dust extinction (rdust.star = 0.6). 

The limiting magnitude constrained by the observa¬ 
tional merger rates basically traces the knee of the LF in 
the K band (see fig. 6 in EMIO) and limits the model pre¬ 
dictions of GNCs to K < 18 mag. The EMIO model cannot 
therefore provide any estimates on the slope change in the 
7f-band GNCs beyond K ~ 18 mag by construction. This 
limiting magnitude is equivalent to the following limits in 
other bands: U 25, B ^ 25, V ~ 25, R ~ 24.5, 7 ^ 24, 
a rsj 23.5, J rsj 20.5, and 77 ~ 19.5 mag. Nevertheless, we will 
show that it complements the simple scenario proposed in 
the EM06 model and provides very interesting clues to the 
nature of the NIR slope change when compared with the 
observational diagnostics derived in the present study. 

According to this model, the bulk of the local massive 
E-SO’s underwent advanced stages of major merger events 
that assembled them at 0.6 < 2 < 1.2, which means that 
these mergers started at 0.8 < 2 < 1.5 approximately. Note 
that the EM06 and EMIO models share two key ingredients: 
1) the massive E-SO’s were definitively assembled in the uni¬ 
verse at late epochs [z < 1.5), and 2) they had experienced 
stages of extreme dust extinction in the early stages of their 
assembly. However, the EM06 model does not assume any 
mechanisms to explain this build-up, whereas the EMIO di¬ 
rectly tests whether major mergers since 2 ~ 1.5 could have 
built them up or not. We remark that the model can fol¬ 
low the evolution of Hubble types later than E-SO’s at all 
masses, but only of E-SO’s which have M* > IO^^Mq at 
2 = 0. 


4 RESULTS AND DISCUSSION 

In i )4.1l we first comment on the results concerning the na¬ 
ture of the 77-band slope change that can be directly derived 
from the observed GNCs by morphological types and of E- 
SO’s by redshift bins obtained in this study. We compare 
these and other GNC data from the literature with the pre¬ 
dictions of different galaxy evolutionary models in 114.21 

4.1 Observational results 

4.1.1 GNCs by galaxy types in the K band 

In Fig.[2l we represent the GNCs in the K band obtained 
in the present study for total, E-SO’s, red highly distorted 
galaxies, red regular discs, and blue galaxies (see Table[T]). 
We have overplotted the GNC compilation of complete data 


performed in the band by CH03, in order to check the relia¬ 
bility of our results by comparing with independent studies. 
Our total GNCs computed for galaxies at 0.3 < z < 1.5 
agree pretty well with the data obtained by other authors 
for magnitudes brighter than our completeness limit (77 ~ 

18.5 mag) down to 77 ~ 16 mag, where our sample starts to 
suffer from poor statistics owing to its area. Therefore, the 
figure indicates that the total GNCs at 16 < 77 < 18.5 mag 
are made up by the galaxy populations at 0.3 < z < 1.5. 

The figure also shows that the galaxy types that ex¬ 
perience a change in the slope of their GNCs around 77 ~ 

17.5 mag are E-SO’s and red discs. However, the contribu¬ 
tion of red discs to the total GNCs is negligible at these 
magnitudes, meaning that the slope change in total 77-band 
GNCs is caused by a change in the trend of the GNCs of 
E-SO galaxies. This decrease in the contribution of E-SO’s 
to the total GNCs at 77 ^ 17.5 mag contrasts with the ris¬ 
ing contributions of blue galaxies and red distorted galaxies, 
which preserve their rising slope up to 77 ~ 18.5 mag. 

Figure[2] demonstrates that the slope change at 77 ~ 

17.5 mag is due to a change of the galaxy population that 
numerically dominates the total 77-band GNCs: from the 
quiescent regular E-SO’s at 77 < 17.5 mag to blue discs at 
77 > 17.5 mag. This is the first direct observational evidence 
for the fact that E-SO’s are the galaxy population respon¬ 
sible for the slope change (as proposed by so many studies, 
see m, because the slope of the total GNCs decreases at 
this magnitude owing to a flattening of the contribution of 
E-SO’s at 77 > 17.5 mag (which slightly decreases, instead 
of continuing to rise, as observed at brighter magnitudes). 

4 . 1.2 GNCs of E-SO’s by redshift bins in the 77 band 

In order to find out the range of redshifts of the E-SO’s 
responsible for the slope change, we have represented the 
77-band GNCs of E-SO’s by redshift bins in Fig.[21 The top 
panel shows the distributions for narrow redshift bins (data 
in Table[2|), while the bottom panel represents the results 
obtained for the three wide overlapping redshift bins (see 
Table[31). We have overplotted the total contribution of E- 
SO’s at 0.3 < 2 < 1.5 from Fig. [2] as a reference in both 
panels. The large error bars in the top panel prevent us 
from deriving any conclusions. However, the GNCs of E- 
SO’s in the wider redshift bins of the bottom panel have more 
reasonable errors that enable us to derive some conclusions. 

The bottom panel of Fig. [3] shows that the slope of 
the total GNCs of E-SO’s at bright magnitudes decreases 
at 77 ~ 17.5 mag because the contribution of these galax¬ 
ies flattens or decreases at 77 > 17 for all redshift bins. 
The E-SO’s that should dominate the GNCs at fainter mag¬ 
nitudes according to the extrapolations of their trends at 
brighter magnitudes are those located at 0.6 < 2 < 1.0 
and 0.8 < 2 < 1.2. However, their contributions decrease 
instead of contributing to raise the GNCs at fainter mag¬ 
nitude s. This is consistent with the results by iBarro et al.l 
ll2009ll . who claimed that the galaxies responsible for the 
slope change were at 2 < 1 (although these authors did not 
identify the type of galaxies). We will show that these GNCs 
of E-SO’s by redshift bins provide evidence of a noticeable 
number evolution of massive E-SO’s at 2 > 0.6 when com¬ 
pared with the models (see il4.2.3ll . 

In conclusion, the 77-band morphological GNCs for dif- 
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Table 2. iC'-band GNCs of red regular spheroid-dominated galaxies (E—SO’s) in narrow redshift bins 



0.3 < 

z < 1.1 

5 

0.4 < 2: < 0.6 

0.6 < z < 0.8 

0.8 < z < 1.0 

1.0 < z < 1.2 

K* 

log Ail 


A* 

logV 

Ai 

Au 

logN 

Ai 

Au 

logN 

Ai 

Au 

log TV 

Ai 

Au 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

(15) 

(16) 

16.25 

2.79 

0.22 

0.22 

2.41 

0.32 

0.32 

1.93 

0.56 

0.52 

2.11 

0.51 

0.43 

- 

- 

- 

16.75 

2.82 

0.21 

0.21 

2.23 

0.41 

0.38 

2.53 

0.32 

0.28 

1.81 

0.83 

0.58 

- 

- 

- 

17.25 

3.15 

0.15 

0.17 

2.63 

0.25 

0.26 

2.71 

0.24 

0.24 

2.11 

0.41 

0.44 

2.33 

0.48 

0.35 

17.75 

3.21 

0.14 

0.16 

2.37 

0.31 

0.34 

2.90 

0.21 

0.21 

2.69 

0.24 

0.24 

1.63 

0.74 

0.70 

18.25 

3.11 

0.16 

0.17 

2.53 

0.28 

0.28 

2.44 

0.33 

0.31 

2.50 

0.26 

0.30 

2.47 

0.36 

0.29 

18.75^ 

3.09 

0.16 

0.17 

1.63 

0.70 

0.72 

2.73 

0.32 

0.23 

2.58 

0.25 

0.27 

2.17 

0.53 

0.41 


* i^-band apparent magnitudes in the Vega system. 

^ The GNCs are provided in units of mag“^ deg“^. Their logarithmic values are listed. 

^ The Ai and Au errors are defined as Ai = log(V) — log(7V — cri) and Au = log(V -h cju) — log(V), where cr^ and (Tu correspond to the 
lower and upper errors due to statistical counting, classification, and cosmic variance uncertainties in each case. 

^ The data in this last magnitude bin are not complete in detection. 


Table 3. iC-band GNCs of red regular spheroid-dominated galaxies (E—SO’s) in wide overlapping redshift bins 



0.3 < 

z < 1.5 

0.4 < z < 0.8 

0.6 < z < 1.0 

0.8 < z < 1.2 

K* 

log Aft 


A* 

logAf 

Ai 

Au 

logAf 

Ai 

Au 

logAf 

Ai 

Au 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

16.25 

2.80 

0.27 

0.20 

2.53 

0.35 

0.26 

2.33 

0.46 

0.32 

2.17 

0.57 

0.37 

16.75 

2.82 

0.26 

0.20 

2.71 

0.29 

0.22 

2.61 

0.32 

0.24 

1.81 

1.67 

0.54 

17.25 

3.15 

0.17 

0.16 

2.97 

0.21 

0.18 

2.80 

0.23 

0.21 

2.50 

0.29 

0.28 

17.75 

3.21 

0.17 

0.15 

3.01 

0.19 

0.18 

3.11 

0.19 

0.16 

2.73 

0.28 

0.22 

18.25 

3.11 

0.19 

0.16 

2.79 

0.26 

0.21 

2.78 

0.25 

0.21 

2.79 

0.26 

0.20 

18.75^ 

3.10 

0.20 

0.16 

2.77 

0.28 

0.20 

2.96 

0.22 

0.18 

2.73 

0.25 

0.22 


* i^-band apparent magnitudes in the Vega system. 

^ The GNCs are provided in units of mag“^ deg~^. Their logarithmic values are listed. 

^ The Ai and Au errors are defined as Ai = log(V) — log(7V — cri) and Au = log(7V -|- ctu) — log(V), where era and (Tu correspond to the 
lower and upper errors due to statistical counting, classification, and cosmic variance uncertainties in each case. 

^ The data in this last magnitude bin are not complete in detection. 


ferent galaxy types and the GNCs for E-SO’s in redshift bins 
evidence that the slope change in the NIR GNCs is caused by 
a change in the contribution of E-SO’s at 0.6 < 2 < 1.2 with 
respect to the expectations derived from the extrapolations 
of their trends at brighter magnitudes. We note that many 
studies had previously proposed this fact (see m, but this 
is the first time that it is shown directly with observational 
data. 


4.2 Clues derived from the comparison with 
evolutionary models 

In this section, we compare the GNC data derived in the 
present and other studies with the predictions of the models 
commented on in We show that a number of conclusions 
can be derived in this way concerning how rapidly the mas¬ 
sive end of the red sequence was established during the last 
-9Gyr. 


4-2.1 Total GNCs from the NUV to the NIR: 
evidence of a relatively late definitive 
build-up of L ~ L* E-SO’s 

Figure|3] compares GNC data in several broad bands from 
the NUV to the NIR together with the predictions of the 
four evolutionary scenarios indicated in fjS] The predictions 
of the PLE, EM06, and EMIO models are plotted for all 
bands. The KW07 model is plotted only for the U, B, R, 
I, and K hlters for clarity. The hgure uses the compilation 
of published GNCs performed by N. Metcalf^^ for the R, I, 
and H bands and for several helds. The data and the equa¬ 
tions used to convert between bands and the AB and Vega 
systems are detailed in his webpage. In the K band, we have 
overplotted the total GNCs obtained in this study, as well 
as the dataset of complete GNCs data compiled by CH03. 
The selection of complete data reduces the typical disper¬ 
sion between authors owing to different extraction proce¬ 
dures, differences in the hlters, and/or cosmic variance be¬ 
tween the helds, which usually smudge the slope change at 

® Galaxy number counts page by the Extragalactic Astron¬ 
omy and Cosmology Research Group of the Durham University: 
http: / / star-www.dur.ac.uk/~nm/pubhtml/counts / counts.html 


© 2000 RAS, MNRAS OOO.nMTI 

















10 


Prieto & Eliche-Moral 



Figure 2. Observational GNCs in the K band obtained in this 
study for all galaxies (total), E—SO’s, blue galaxies, red highly 
distorted galaxies, and red regular discs (see the legend in the 
figure). Vertical dotted line', limiting magnitude for 100% com¬ 
pleteness of our data. Black filled circles: complete GNCs dataset 
compiled by CH03 in the K band from the literature, plot¬ 
ted for compa r ison llGardner et alj|l993l: IPiorgovski et al.lll9^: 
McLeod_et_ajj_[ljl95ji__[^ian^^r^aribe^ Jl997l : iMinezakieT alj 
igO^ T ISzokoIvet al.lll99'ST^tanr'er'^l200ll . CH031. Solid lines: 
fits performed to the bright and faint ends of the GNCs by CH03. 
The slopes obtained by these authors are indicated in the figure. 
Black arrow : location of the slope change in total GNCs. [A colour 
version is available at the online versioni\ 


K ~ 17.5 mag. We have also re presented the deep GNCs ob¬ 
tained bv iRetzlafF et al] (120101 ) in UDF and GOODS-Sonth 
fields in the J, H, and K bands. References for all data in 
the figure are provided in the caption. 

The PLE model cannot fit the GNC data at faint mag¬ 
nitudes in any band, as shown in Fig. 3] The predictions of 
PLE are above the data in the K band and below them 
in blue optical filters. In particular, the PLE model rolls 
over at R' ~ 20 magnitude, basically because of the evolu¬ 
tion of the cosmological volume element, the change in the 
k-corrections in the K band from being positive to negative, 
and the different shape s of the LFs representing E-SO’s and 
spirals (lBershadvll2003l) . However, this slope change appears 
to be ~2.5mag deeper than the observational slope change 
at if ~ 17.5 mag. This mismatch indicates that the PLE 
model predicts too few blue galaxies and too many red ones 
at high redshifts, which clearly indicates the existence of 
a conversion of blue into red galaxies in the universe (see 
references of previous studies already stating this in ©• 

The KW07 model fits the observations pretty well in 
the U, B, R, and I bands, but clearly fails to reproduce to¬ 
tal GNCs in the K band faintwards of if ~ 19 mag. This 
discrepancy was already noticed by these authors, who con¬ 
cluded that the evolutionary scenario of the Millenium sim- 




Figure 3. Observational if-band GNGs of E-SO galaxies ob¬ 
tained in the present study for different redshift bins up to z < 1.2 
(see the legend in the figure). The total GNGs of E—SO’s in the 
band for 0.3 < z < 1.5 from Fig.[2]have been overplotted for refer¬ 
ence. Top panel: for narrow redshift bins. Bottom panel: for wide, 
overlapping redshift bins. Vertical dotted line: limiting magnitude 
for 100% completeness of our data. Black arrow: location of the 
slope change in the total GNGs. [A colour version is available in 
the online version] 

ulation ‘overpredicts the abundance of moderately massive 
galaxies at high redshift, despite the fact that late merg¬ 
ing plays a major role in the build-up of its more mas¬ 
sive galaxies’. This statement means that massive galax¬ 
ies are being assembled at earlier times in the model (in 
fact, at z > 2) than in the real universe. This model may 
be missing some physical mechanisms that delay the defini¬ 
tive major assembly of these galaxies until later epochs (at 
z < 1. 5), according t o the total GNCs from NUV to NIR 
bands. iGrazian et al.l (l2009l ) also compared deep U and K 
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Figure 4. Predictions of the models on the total GNCs from the NUV to the NIR bands compared with real data {coloured circles). 
The GNCs of each band has been displaced in the magnitude axis by a constant for clarity, as shown in the legend of data colours at 
the top. Black solid lines: predictions of the EMIO model in all bands. Green dotted-dashed lines: predictions of the EM06 model in all 
bands. Light-blue long-dashed lines: predictions of the KW07 model in the f/, R, R, 7, and K bands. Grey dashed lines: predictions of 
the PLE model for all bands. This figure uses the compilation of published GNC data by N. Metcalfe in the R, 7, and 77 bands based 
on CCD imaging. The GNCs for the K band include our total GNCs, the selection of complete data performed by CH03, and the deep 
data bv iRetzlaff et ahl l|2Ql(l) . 


The r eferences for each band are the following: U ban d <lSorLg;mla_et__^ [l9^; [AriiQuts_et__^ [200l|; h^jSud^^^^ 


Capak et alJ 


dArnouts et al.l 120011: 

Huaner et al.l 120011: iKiimme 

& Waenet! I 2 OOII; 

Yasuda et al.l 120011: iMcCracken et al.l 1200251: ICapak 

et al.1 12004 

iKashikawa et al.1 

2004 

: lEliche-Moral et al.1 l2006bl: 

Rovilos et al.112009 

). V banddCardner et alJ Il996l: 

Capak et al 

|2004 

), R band 

dCouch & Newell! 

1984 

:lHall &; Mackav 

ll984l:IStevensonlll985l:llnfante et al.lll986l: lKodll986l: lYee &: Green! 

19871: iTvson 

198J 

Jones et al 


1991 

; 

1995 

; 

2001 


1986 


1995 

; 


; Bertin &; Dennefeld Il997l: iHogg et al.lll997l:lMetca?fee^an^^8l:lArnoutsetS.lll99^.l200lljHuan^^^an[20^|jK(hninel^^h!gn^ 


2004), 7 band <|Hah fc^ac^^ 


Metcalf^talJ 200l|j McCrad^en et aI] [200j; Ca£al^^ayj2004 


■ .Kashikawa et al.l 

I1995I: IPriver et ^ 


Koc 


199^ jGiazebrooket^^^^^|j^e^]e^eet^L 


Carder et alj|l996l: Arnou^ et al.lll999l . |20^ j_ Metcalf^et*aLn200ll: lYasud^T^ahf 

'Capa^et alJ I2OO4I . I2OO7I : iKashik^^et a^ l2 Q04h. .s band <|Capak ^^^i?TT20Q4h. J bai^ llBershadvet^ll 


Smail et al 


McCracken et al.l 


199^ 


2002 


199d: IVhishn^ et alj 


2000LJS^^^^^^^ay[|20^|^^ni^^^J' 200^ ICristQbal-HornillQr^^^.l 1200^: 


Retzlaff et al 


l|Teplitz^^ar|~[l99^: |Yan et al.l |l99^: |Thompson et al.| |l999|: |Ma rtini | I 200 J: iMetc^e et al. | |2001 : I Chen e^aLTI^OOS 

[Cristbba^Homillos et a?Tl2Q0^: lKeena^^^anr2010l: iRetzlaff e^alTI^^^ . K band 


_■ _ _ __ _ _ __.Frith et al 1 l2006l: 

iRetzlaffe^^ I 2 OI 0 I . results derived in the present 


201 


K 


H band 


study, and references in the caption of Fig.[ 2 jl- [A colour version is available in the online version.] 
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GNCs with the predictions of three hierarchical CDM mod¬ 
els (including also the KW07 model, see their figs. 3 and 4). 
They found large discrepancies at the faint end of the U- 
band GNCs between the three models that clearly favour 
the KW07 model. They attributed these discrepancies to 
differences between the models in the recipes dealing with 
the star formation activity and dust extinction of the faint 
galaxy population a.t z < 2. 

The EM06 and EMIO models fit well all the data 
from NUV to NIR in their ranges of validity. This good¬ 
ness of fit suggests that they may be describing the build¬ 
up of the bulk of the E-SO’s more realistically than the 
PLE and KW07 models in global terms. In particular, 
the EM06 model, besides reproducing the slope change at 
K 17.5 mag, can also predict the faint end of the GNCs 
in the K band up to K ~ 24 mag drawn by the data 
by iRetzlaff et al] (l2010ll . These authors reported that they 
found no change in the slope of the K GNCs at this very 
faint end (i.e. for K > 20mag), so it is extremely significant 
that the EM06 model reproduces it, whereas the PLE model 
cannot. Both the PLE and EM06 models assume PLE for 
all Hubble types later than Sa according to standard SFHs 
and very high formation redshifts {zf = 3), so, in this sense, 
they are analogous. They differ basically in that the PLE 
model assumes that E-SO’s appear on the cosmic scene at 
Zf = 3, whereas the EM06 model delays their appearance 
until Zf = 1.5. Therefore, the shape of the total Tf-band 
GNCs directly rejects passive evolutionary scenarios for E- 
SO’s from very high redshifts not just because of the slope 
change a.t K ^ 17.5 mag, but also because of their feature¬ 
less slope at the very faint end. 

Although the EMIO has a much shorter magnitude 
range of validity (see Em, it nearly follows the predictions 
of the EM06 model, in particular, the slope change until 
K — 18 mag. We note that the EM06 and EMIO models dif¬ 
fer from the PLE and KW07 models in that they propose a 
late build-up of the bulk of the massive early type galaxies 
(at 2 < 1.5). The EM06 model imposed it ad hoc, whereas 
the EMIO model derived it naturally from the major mergers 
registered up to 2 ~ 1.5 (assuming that one E-SO has risen 
from each observed major merger). In this sense, the EM06 
model can be considered a coarser approximation to the sce¬ 
nario proposed in the EMIO model, but with the clear advan¬ 
tage that it can make predictions at all magnitudes, whereas 
the EMIO model is limited to Mb < —20 mag. Therefore, 
we can complement the predictions of the EMIO model by 
those of the EM06 model at very faint magnitudes, consider¬ 
ing the latter to be a rough sketch of the formation scenario 
proposed in the former. 

Note also that the EM06 and EMIO models are capable 
of reproducing ?7-band GNGs down to ?7 ~ 26-27 mag at 
least as well as the KW07 model, whereas the PLE model 
fails to predict them at the faint end. Therefore, models 
assuming PLE evolution of E-SO’s since very high redshifts 
can be neglected not only by attending to the GNCs in the 
K band, but also to the deep GNCs in the optical-to-NUV 
bands. 

Figure[5] shows the contribution of the various galaxy 
types to the total GNCs predicted by the EMIO model in 
the U, B, and K bands, compared with the same observa¬ 
tional data as in Fig.[4l The model is capable of reproducing 
the slope change until A ~ 18 mag (the limiting magnitude 



CVi 



Figure 5. EMIO model predictions of the contribution of the 
various morphological types to the total GNCs in the f/, and 
K bands compared to observational data (see references in the 
caption of Fig.[4|l. The region where the predictions of the EMIO 
model are not valid in each band is shaded in light blue. Black 
thick solid lines: total GNCs predicted by the EMIO model in each 
band. Grey dashed lines: predictions of the PLE model on the to¬ 
tal GNCs in each band, plotted as a reference. Rest of lines: EMIO 
model predictions of the GNCs in each band due to the E-SO’s 
(red solid), Sa-Sab’s (yellow dashed), Sb-Sbc’s (olive green dotted- 
dashed), Sc-Scd’s (dark green three dotted-dashed), Sd-Irr’s (blue 
dotted), and red advanced stages of major mergers (purple thick 
long-dashed). [A colour version is available at the online version.] 


of the model) and attributes it to a flattening of the con¬ 
tribution of E-SO’s to the total GNCs at K > 16 mag. The 
GNCs of the later types (from Sa’s to Sd’s) do not signifi¬ 
cantly change the rising trend they exhibit at bright magni¬ 
tudes after the break. The figure shows that i^-band GNCs 
are controlled by E-SO’s up to the slop e change. Accounting 
for the shape of their LF up to 2 : ^ 2 (|Arnouts et al.ll2007l : 
ICirasuolo et al.ll2007li . NIR GNCs are controlled by L > L* 
E-SO’s up to A ~ 17.5 mag according to this model. We 
will also show that the addition of the GNCs from spiral 
types in the model exceeds the GNCs of E-SO’s after the 
break, so the global scenario proposed by the EMIO model 
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is in good agreement with the observational results that we 
have derived in l|4.1l on the nature of the slope change in the 
ff-band GNCs (see 34101) . 

To summarize, the observed distributions of total GNCs 
from the NUV to the NIR bands support a late definitive 
build-up for the majority of the E-SO population (at z < 1.5) 
against the scenarios proposing that they have evolved pas¬ 
sively since very early epochs (at z > 2). This late epoch of 
assembly of massive E-SO’s (postulated by the EM06 model 
and predicted by the EMIO model under the assumption 
that each observed major merger up to z ~ 1.5 has produced 
an E-SO galaxy) is in good agreement with traditional hier¬ 
archical scenarios of galaxy form ation llPe Lucia et al ]|2006l. 
120071: 157 Lucia fc Blaizot lboOTh. Bu t it also agrees with the 
recent results bv IChoi et al.l (1201411 . who report that the 
light-weighted ages of massive galaxies at 2 < 0.7 indicate 
that they are noticeably younger than the age of the uni¬ 
verse at all redshifts up to 2 ~ 0.7, which would suggest 
an effective single-burst star-formation epoch of 2 < 1.5 for 
them. 


4-2.2 GNCs by redshift bin: evidence of the 

progressive assembly of massive E—SO’s at 

0.8 < 2 < 1.5 

Figure[6] represents the predictions of the EM06 and EMIO 
models on the ff-ban d GNCs by r edshif t bins, compared 
to the data derived bv iBarro et al.l ll2009l L The predictions 
of the PLE model are shown in both cases for reference 
because the PLE scenario for the evolution of E-SO’s since 
high redshifts is rejected when attending only to the total 
GNCs from NUV to NIR bands In fact, the PLE 

model adequately fits real data up to 2 ~ 0.75 in the figure, 
but it lies above the observations at 2 > 0.75, and predicts 
too many luminous galaxies and too few faint ones at 1.25 < 
2 < 2. This model lacks a physical mechanism that converts 
massive galaxies into several low-mass components as we 
counter-evolve galaxy populations backwards in time. 

The predictions of the EM06 model (top panels of the 
figure) are in better agreement with observations for L > L* 
galaxies than the PLE model, but they fail to reproduce the 
faint end of the GNC distribution from 0.5 < 2 < 0.75 
back to earlier epochs. Therefore, assuming that most E- 
SO’s formed at 2 ~ 1.5 provides a good approximation to 
the evolution of the bright end of the GNCs by redshift bin, 
but it does not work for the faint end. This suggests that the 
build-up of this galaxy population must have taken place at 
different epochs at 2 < 1.5, i.e. that the assembly of these 
galaxies at 2 < 1.5 has been progressive. 

In the middle panels of Fig.[6]we show that, within its 
limitations in magnitude range, the EMIO model reproduces 
the evolution of the massive end of the A'-band GNCs by 
redshift bin. The growth of massive E-SO’s is progressive 
in this model and basically takes place at 0.8 < 2 < 1.5, 
whereas in the EM06 model (top panels) all these galaxies 
are inserted into the cosmic scenario at 2 = 1.5. However, the 
progressive build-up of the massive E-SO’s at 0.8 < 2 < 1.5 
in the EMIO model describes the bright end of the total 
GNCs by redshift bin in the K band in a similar manner as 
the sudden build-up of these galaxies in the EM06 model. 
We explain the reason for this below. 

We have also overplotted the contribution of the vari¬ 


ous galaxy types to the total GNCs expected by the EMIO 
model at each redshift bin in the middle panels of Fig.|6] The 
bright end of TF-band GNCs is controlled by the E-SO’s up 
to 2 ~ 0.75, but another population (those galaxies under¬ 
going advanced stages of major mergers) starts to contribute 
similarly to the total GNCs from this epoch up to 2 ~ 2 in 
the model. Note that the addition of the GNCs of E-SO’s 
and these advanced stages of major mergers at 0.75 < 2 < 2 
controls the bright end of the GNCs by redshift bin in the K 
band. The difference from the EM06 model in the top pan¬ 
els is that the EM06 model does not consider the existence 
of galaxies undergoing major mergers. However, the total 
Lf-band GNCs of E-SO’s and galaxies in advanced stages of 
major mergers predicted by the EMIO almost coincides by 
chance with those of E-SO’s expected by the EM06 model 
at all redshift bins, and they both reproduce the observed 
data similarly. Therefore, the EM06 model attributes the 
evolution observed in the bright end of the 7F-band GNCs 
by redshift bin as a mere colour evolution of the E-SO’s up 
to 2 1.5 and to their vanishing at earlier epochs, whereas 

the EMIO model explains it through a progressive build-up 
of the massive E-SO’s at 0.8 < 2 < 1.5 through the major 
mergers reported strictly by observation. 

In the bottom panels of Fig.[^ we show the difference 
be tween the predictio ns of the three models and the data 
bv iBarro et al.l (l2009l f in dex as a function of the K magni¬ 
tude. The numbers in each panel correspond to the ratios of 
the total number of galaxies predicted by each model with 
respect to the total number of galaxies according to the ob¬ 
servational data, integrated down to the limiting magnitude 
set by the EMIO model in each redshift bin. Attending to 
the distribution of differences in dex, the predictions of the 
EM06 seem to fit the data better than the EMIO and PLE 
models in the range of validity of the EMIO model, but when 
the integrated ratios are compared, the EMIO model pro¬ 
vides better predictions globally than the EM06 model at 
2 > 1. The EM06 model also predicts Lf-band GNCs higher 
than those expected by the PLE model by 0.3-0.5 dex at 
faint magnitudes for 2 > 1. 

The results of the EMIO model agree pretty well with 
observational evidence that massive galaxies have experi¬ 
enced passive evolution o nly since 2 ^ 0 .7 according to 
their abundances of metals jChoi et al.l2014li . These authors 
claim that newly quenched galaxies are added at 2 < 0.7 
only at the lower masses (log(M,/MQ) < 10.5), which 
means that the bulk of the build-up of the more massive 
ones must have taken place at 2 > 0.7 (in agreement with 
the predictions of the EMIO model). This growth of the qui¬ 
escent populatio n only at the lowest masses is also supported 
by other studies dHahn et al.ll2014^ . However, other investi¬ 
gations find evidence of that the brightest quiescent mem¬ 
bers continue to acc rete mass throu gh dry major mergers in 
clusters at 2 < 0.6 (iLiu et al.ll2015l l , so the real scenario of 
the build-up of these galaxies must be more complex. 

In summary, Lf-band GNCs by redshift bin are an ob¬ 
servable that directly discards the PLE model at 2 > 0.8 
for massive E-SO galaxies and the ‘sudden’ appearance 
of the whole population at 2 ~ 1.5, as proposed by the 
EM06 model (at least at the faint end). Instead, these GNC 
data imply a progressive build-up of massive galaxies at 
0.8 < 2 < 1.5. 
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Figure 6. Predict ions of the models of the total K-h&nd GNCs in different redshift bins up to 2 : = 2, compared to the real data from 
iBarro et alj ll2009h . Top panels: data compared to the PLE and EM06 models. Middle panels: data compared to the PLE and EMIO 
models. We also represent the contributions to the total GNCs in each redshift bin by morphological types expected by the EMIO model. 
The legend for these lines is the same as in the bottom panels of Fig. (5] Bottom panels: differences in dex between the total GNCs in 
the K band predicted by the PLE, EM06, and EMIO models and the observational data in each redshift bin as a function of the apparent 
iC-band magnitude. The factors written in each panel at the bottom represent the ratios of the total number of galaxies predicted by 
each model (PLE, EMIO, and EM06) with respect to the total number of galaxies according to the observational data, a ll integrated 
down to the limiting magnitude set by the EMIO model at each redshift bin. Yellow squares: observational data obtained bv iBarro et alj 
1 I 2 OO 9 I L The typical observational errors at three magnitudes spanning the observational range have been overplotted, but displaced 
vertically from their data points for clarity. Grey dashed lines: PLE model. Green dotted-dashed lines: EM06 model. Black solid lines: 
EMIO model. Region shaded in light blue: Magnitudes for which the predictions of the EMIO model are not valid for the total GNCs in 
the K band. [A colour version is available in the online version. 


4-2.3 GNCs of E—SO’s by redshift bin: this 

assembly involved 50% of present E—SO’s at 
least 

In Fig. [7] we compare the GNCs of E-SO’s by wide redshift 
bins obtained in 114.1.21 with the predictions of the EMIO 
model of the GNCs of E-SO’s. The expectations for PLE 
are shown for reference. We also show the percentages rep¬ 
resented by the total number of galaxies predicted by the 
EMIO model and the data with respect to the expectations 
of the PLE model for the magnitude range in which the 
EMIO model is valid. This hgure shows that the EMIO model 
provides predictions much more in agreement with the data 


than PLE. The GNCs of E-SO’s predicted by the EMIO 
model at the bright end are a factor of ~2-3 times lower 
than the prediction of the PLE model at 0.8 < 2 ; < 1.2. 
Comparison of the real data with the PLE model implies 
that the number density of massive E-SO’s is lower than 
their local value by a factor of ~2 at 0.6 < 2 < 1 and 
at 0.8 < 2 < 1.2. Similar values are predicted by the 
EMIO model. The GNCs of E-SO’s by redshift bin reveal 
that more than ~50% of present-day massive E-SO’s have 
been progressively built up at 0.8 < 2 < 1.5. The EMIO 
model attributes most of this assembly to the major merg- 
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ers reported by observations (mainly through wet ones at 
high redshifts, see EMIO). 

The Tf-band GNCs of E-SO’s by redshift bin thus sup¬ 
port the key prediction of the EMIO model that the pop¬ 
ulation of massive E-SO’s decreases down to -^40-50% at 
0.6 < z < I with respect to the PLE extrapolation of their 
local population and down to 20- 30% at 0.8 < a < 

1.2. This numerical d ecrease of massive E-SO’s agrees well 
with the estimates b y ICristobal-Hornillos et al.l (l2009l l and 
iMoresco et al.l (120131 '). This finding implies that we may con¬ 
sider the bulk of the present-day massive E-SO population 
to have been definitively assembled at z ~ 0.8-1, so that 
they may be considered to have been in place since 2 ~ 0.6 
(i.e. only during the last 6 Gy r of cosmic histo ry), in 
good agreement with the results bv IChoi et al.l (l2014l ) noted 
above. 

The EMIO model attributes this build-up entirely to 
major mergers. If this hypothesis turned out to be approxi¬ 
mately true, it would mean that massive galaxies have expe¬ 
rienced a growth in stellar mass of a factor x2 since 2 ~1- 

1 .2, which is in good agreement with the recent estimates by 
IShankar et al.l ( 20151 ). These authors also conclude that hi¬ 
erarchical models could explain this growth through major 
mergers, and that the rate expected for them is very similar 
to observational estimates, thereby supporting the scenario 
proposed by the EMIO model as well. 

In summary, the /S'-band GNCs of E-SO’s by redshift 
bin support the claim that present-day massive E-SO’s expe¬ 
rienced a noticeably progressive assembly at 0.8 < 2 < 1.5, 
involving at least ~50% of their present-day number den¬ 
sity. The EMIO model shows that it is feasible to explain 
this build-up just through the major mergers reported by 
observations up to 2 ~ 1.5. 


4-2.4 GNCs by morphological type in different 

bands: major mergers as possible drivers of 
this assembly 


In Fig.|5] we show the GNCs in the I band by morphological 
type obtained by various authors, compared to the predic¬ 
tions of several models. We have overplotted the expecta¬ 
tions of the PLE and EMIO models, as well as the GNC 
predictions by galaxy type deriv ed bvlBaueh et al.l l|l996bl ) 
from their semianalytic model jBaugh et al. 1996ah . This 
model predicts the galaxy evolution starting from the power 
spectrum of primordial density fluctuations, and implement¬ 
ing the effects of different physical processes ruling galaxy 
formation and evolution, such as gravitational collapse, gas 
cooling, star formation, feedback, and galaxy mergers. This 
model adequately reproduces the slope of the t otal GNCs in 
the I b and, but not in the B band (see Fig. 1 in iBaugh et al.l 
Il996tj) . 


Figure|8] demonstrates that GNCs by morphological 
type in optical bands are not very sensitive to galaxy evolu¬ 
tion. Independently of the data dispersion in the hgure, the 
three overplotted models (which are very different among 
themselves) provide similar predictions for the GNCs of E- 
SO’s, spirals, and Irr-|-Merger galaxies in the I band. There¬ 
fore, little information can be obtained from morphological 
GNCs in optical bands. 

However, we have found that morphological GNCs in 
the NIR bands are much more sensitive to the galaxy evo¬ 


lution experienced by each galaxy type than those in the 
optical bands. Figure|9] shows the total GNCs by morpho¬ 
logical types in the K band derived in the present study 
f M.l.lll compared to the predictions of the EM06 model 
(left panel) and the EMIO model (right panel). We have 
had to establish an ad hoc analogy between the data types 
defined in P13 (see 31} and those defined in each model. The 
broad galaxy classes assumed by the EM06 model makes it 
difficult to compare the predictions of this model with the 
available data. We have identified the SO/a-Sb type in the 
model with red discs, but it is obvious that this class con¬ 
tains galaxies that are not strictly ‘red’ (the Sb’s), which 
should be included in the comparison with real blue galax¬ 
ies instead. This explains why the model overpredicts the 
GNCs of red discs and underestimates the GNCs of blue 
galaxies at all magnitudes. 

We have shown that the simple scenario of the EM06 
model can be rejected simply by attending to its predictions 
at the faint end of the total if-band GNCs by redshift bin 
(see Fig. [6}. The disagreement between the data and the 
model predictions was interpreted as a sign that the bulk of 
low-mass {L < L*) E-SO’s cannot have formed at a short 
time period around 2 ~1.5. Their definitive assembly was 
probably more gradual over the last ~9 Gyr. Attending now 
to the predictions of this model for the GNCs of E-SO’s in 
Fig.m we can also refute a sudden formation of the bulk 
of bright E-SO’s (L > L*) at 2 ~ 1.5. The model clearly 
overpredicts their contribution at if > 16.5 mag. Moreover, 
the EM06 model does not provide any explanation for the 
existence of red highly distorted galaxies found in real data. 
It postulates the existence of initial stages of extreme dust 
extinction in the E-SO’s after their appearance on the cosmic 
scene. But these galaxies remain E-SO’s at those stages, they 
do not necessarily exhibit distorted morphology. 

Nevertheless, the GNCs by morphological type pre¬ 
dicted by the EMIO model in the K band (right panel of 
Fig.E} can approximately reproduce the observed distribu¬ 
tions of all galaxy types identified in this study up to the 
limiting magnitude of validity of the model. To compare data 
and model predictions, we have identihed the observed red 
regular spheroid-dominated galaxies with the E-SO’s in the 
model, real blue galaxies with modelled Sb-Irr’s, red regu¬ 
lar disc-dominated galaxies with Sa-Sab’s, and hnally red 
highly distorted galaxies with the advanced stages of major 
mergers predicted by the model. The EMIO model is capable 
of explaining the shape and numbers of the GNCs of mas¬ 
sive E-SO’s up to if ~ 18 mag, as well as the rise in GNCs 
of red distorted objects at faint magnitudes. Contrary to 
the EM06 model, the EMIO model naturally predicts high 
levels of dust extinction in these early phases of the assem¬ 
bled E-SO’s because they are assumed to derive mostly from 
gas-rich major mergers, which are known to undergo intense 
starbursts at advanced stages of the interaction (see refer¬ 
ences in EMIO). The success of the EMIO model is thus two¬ 
fold: first, it predicts the GNCs of bright E-SO’s observed 
since 2 ~ 1.5, and, secondly, it offers an explanation of the 
appearance of a population of highly distorted galaxies on 
the red sequence at 2 > 0.8 (marginal at lower redshifts), 
as the advanced stages of the major mergers reported by 
observations up to 2 ~ 1.5. This shows that it is feasible 
to explain the assembly of massive E-SO’s at 0.8 < 2 < 1.5 
via major mergers. Note that it can also nearly explain the 
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Figure 7. Predictions of the EMIO model for the GNCs of E-SO’s in wide, overlapping redshift bins in the K band, compared to the 
data derived in the present study (see Fig. mi. Circles: GNCs of E-SO’s for the redshift bin indicated in each panel. Red solid line: GNCs 
of E-SO’s predicted by the EMIO model in the corresponding redshift bin. Dark red dashed lines: GNCs of E—SO’s predicted by the PLE 
model in each redshift bin. Blue shaded region: Region where the predictions of the EMIO model are not valid for total R'-band GNCs. 
Vertical dotted line: limiting magnitude for 100% completeness of our data. The percentages written in each panel represent the fraction 
of the total numbers of E—SO’s derived from the observational data and expected by the EMIO model with respect to the number in the 
case of PLE, integrated from the lowest magnitude sampled by the observational data down to the limiting magnitude set by the EMIO 
model at each redshift bin. [A colour version is available in the online version.] 
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Fi gure 8. Pred i ctions of the PLE an d EMIO models f or the GNCs by morphol ogical type in the I ban d, compared to the data obtained 
bv IPriver et alj ll 19981. green c^rc/es). IPhillipps et a D ll2000l. orange stars), and I Volonteri et al.l ll2000l . yellow triangles). Left panel: for 
the E-SO’s. Middle panel: for the Sa-Sc’s. Right panel: for the Sd—Irr’s and mergers. Solid lines: predictions of EMIO model for the 
galaxy type corresponding to each panel. Dashed lines: predictions of the PL E model corresponding to each type, shown for reference. 
Dotted lines: predictions of the semi-analytical model bv iBaugh et al.l lll996bl ) for each galaxy type. [A colour version is available in the 
online version.) 


GNCs of red disc-dominated galaxies if they are identified 
with Sa-Sab’s. 

The EMIO scenario is consistent with the results derived 
from the observed GNCs of Extremely Red Objects (EROs). 
EROs are considered as the counterparts of local L > 2- 
3 X L* galaxies at 0.7 < z < 1.5, with the bulk of population 


at 2 ~ 1 llVaisanen fc Johanssonir2004l l. In the framework of 
the model, the L > L* E-SO’s derive from blue starform¬ 
ing galaxies through the major mergers that are observed to 
occur at 0.8 < a < 1.5, meaning that these galaxies are un¬ 
dergoing advanced stages of these mergers at 0.6 < z < 1.2, 
probably suffering intense dust extinction that reddens them 
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significantly. Moreover, many other studies support that 
there must have been a significant assembly of red galaxies 
at t he expense of the blue ones during th e last 10 Gyr (see, 
e.g., iFaber et al.ll2007l : iKriek et al.ll2008|. P13~). In particu - 


lar, the model predictions agree with lBuitraeo et al 


. p articu - 


who found that the early-type galaxies have been the pre¬ 
dominant morphological class for log(M«/M 0 ) > 11 since 
only 2 ~ 1. 

The rise of spheroid-dominated galaxies at the expense 
of these highly distorted objects at the massive end of the red 
sequence at 0.8 < z < 1.5 can be observed in the evolution 
of the relative fraction of both types in colour-magnitude 
diagrams. In Fig. [10] we plot the apparent B — K vs. K di¬ 
agrams for our blue and red samples in the three wide red- 
shift bins defined in P13. The various morphological types of 
red galaxies are identified by distinct symbols. The expected 
traces for different galaxy types within each redshift interval 
have been overplotted. They correspond to the same mod¬ 
els plotted in Fig.[T] but different stellar masses have been 
assumed to scale the models in each redshift bin in order to 
overlap them with the data: M* = lx, 2x, and 4x IO^^Mq 
for 0.3 < 2 < 0.7, 0.7 < z < 1.1, and 1.1 < a < 1.5 respec¬ 
tively. As times goes by, the relative presence of distorted ob¬ 
jects in the red sequence decreases, whereas the relative frac¬ 
tion of regular spheroid- and disc-dominated galaxies rises. 
This figure also explains qualitatively why the appearance 
of the slope change in GNCs at A” ~ 17.5 mag. The GNGs 
at A < 17 mag are numerically controlled mostly by the red 
objects at 0.3 < a < 1.1 (compare the relative numbers of 
red and blue objects at A < 17 mag in the three redshift 
bins). However, the number of blue objects at 0.3 < a < 1.1 
tends clearly to dominate them at A > 18.5. 

The EMIO model shows that it is feasible to attribute 
the build-up of massive E-SO’s entirely to major merg¬ 
ers, but this scenario is an obvious oversimplification of 
reality. It does not account for additional evolutionary 
mechanisms, such as gas stripping and st rangulation in 
clusters, simple fading, or harassment (see lAguerril I 2 OI 2 I : 
iKormendv fc Bendeill20l3 . for a review on the possible pro¬ 
cesses), that are known to have played a relevant role in the 
migration of galaxies from the blue cloud to the red sequence 
(mostly producing SO’s). However, the evidence of the pro- 
tagonism of major and minor merging in the formation of 
massive lenticulars has been growing in recent years, so the 
model could be providing a rough (but almost realistic) pre¬ 


ulation fsee. e.E.. Eliche-Moral et al.ll2006al: iBernardi et al. 

2011 b 

lal: lEliche-Moral et al.l201ll. 20121.12013 

: Peirani et al. 

20091: 

Yane et al. 2009 

Hammer et al.ll2009a 

. I 2 OI 2 I. 2009bl: 

Bernard! et al.l 2011cJ 

3 : Borlaff et al.l 20141: Tania et al. 

2014 Choi et al.l 2014 

Ouereieta et al.ll2015l') 

. The scenario 


proposed in the EMIO model is also supported by the recent 
results by I Wagner et al.] (I 2 OI 5 I I. who indicate that the num¬ 
ber evolution of the majority of massive early-type galaxies 
at 2 >1.2 seems to be driv en by recent gas-rich mergers (see 
also IWhitaker et al.ll2012ll . 

In summary, GNGs by morphological type are much 
more sensitive to the evolution of massive E-SO’s when de¬ 
rived in NIR bands than in the optical bands. The morpho¬ 
logical GNGs of E-SO’s and highly distorted objects on the 
red sequence are compatible with a scenario in which the 
build-up of massive E-SO’s at 0.8 < 2 < 1.5 has occurred 


through the major mergers reported by observations up to 
2 ~ 1.5 (in good agreement with hierarchical theories of 
galaxy formation). The GNGs of E-SO’s (total and by red¬ 
shift bins) derived in the present study indicate that the 
slope change in the total GNCs at A ~ 17.5 mag is nothing 
other than a vestige of the assembly of more than a half of 
present-day massive E-SO’s at 0.8 < 2 < 1.5, possibly driven 
by major mergers in their majority. 


5 SUMMARY AND CONCLUSIONS 

We have studied the galaxy population and evolutionary 
processes responsible for creating the slope change registered 
at A ~ 17.5 mag in NIR GNCs. We have derived the con¬ 
tribution of various galaxy types to the total A-band GNCs 
at 0.3 < 2 < 1.5, as well as the GNCs of E-SO’s by redshift 
bin, in order to identify the populations that produce this 
feature. We have used the data and galaxy samples selected 
in P13. The GNCs and other data obtained from the litera¬ 
ture have been compared with the expectations of a number 
of galaxy evolutionary models. 

We show observationally that the slope break in the A- 
band total GNCs is caused by a sharp change at 2 < 1.5 in 
the galaxy population that controls them in terms of their 
numbers (from the quiescent E-SO’s at A < 17.5 mag to the 
blue star-forming discs at fainter magnitudes). The change 
is generated by a substantial flattening of the contribution of 
E-SO’s to the total A-band GNCs at 2 > 0.6, as compared to 
their rising trend at the bright end of the GNCs distribution. 

From the comparison of the GNC data obtained here 
and from the literature with the models analysed, we have 
found the following results: 

(i) We corroborate that models in which massive E-SO’s 
evolve passively since high redshifts (2 > 2 ) cannot predict 
the slope change in the A band GNCs and are discarded by 
attending only to the total GNCs from the NUV to the NIR 
bands. On the contrary, we show that the models (such as 
EM06 and EMIO) that assume a significant assembly of this 
galaxy population at 2 < 1.5 can simultaneously reproduce 
the total GNCs from the NUV to the NIR bands, including 
the slope break at A 17.5 mag. 

(ii) The total A-band GNCs by redshift bin reveal that 
massive galaxies (mostly E-SO’s) have experienced a signif¬ 
icant progressive assembly at 0.8 < 2 < 1.5. 

(iii) The GNCs of E-SO’s by redshift bin demonstrate 
that this assembly at 0.8 < 2 < 1.5 involved at least ~50% 
of the present-day massive E-SO population (L > L*). 

(iv) GNCs by morphological type are more sensitive to 
galaxy evolution when derived in the NIR than in the red 
optical bands. 

(v) The EMIO model, which considers massive E-SO’s to 
have derived from the major mergers registered in obser¬ 
vations, is capable of simultaneously reproducing the total 
GNCs from the NUV to the NIR bands, as well as total 
GNCs of E-SO’s by redshift bin in the A band. It can also 
reproduce the GNCs by morphological type if the following 
identification between the modelled and real types (as de¬ 
fined in P13) is established: the modelled E-SO’s correspond 
to the observed red regular spheroid-dominated galaxies, 
Sa-Sab’s in the model are the regular discs detected on the 


© 2000 RAS, MNRAS 000. OMR 


















































18 Prieto & Eliche-Moral 



be 

<u 


bU 

a 

s 


bo 

O 


5.0 


4.5 


4.0 - 


3.5 - 


3.0 


2.5 - 


2.0 - 


1.5 


..... 

■-PLE - Total 

. 

• 

FMIO - Totnl 


^ • 

EMIO - E- 

-so 


^ 

EMIO - Sa-Sab 


■ _ _ EMIO - Adv. mai. mercers 


_ EMIO - Sb-Irr ^ 


7 Data (present study): 

m 

□ Total 




o Red bulge-dom. 


■ 

O Red disc—dom. 



☆ Red hiffhlv— 



■dr \ 

distorted 


I' 7 , 


A Blue glxs. 


$ : 

X < 

> 





I f 


A y 







" 

/a 

< 



/ /•' < 

> i 

r 


- / A 

i 


- 

r 

/ 








14 15 16 17 18 19 

K (mag) 


20 


21 


Figure 9. X-band GNCs by morphological type derived in the present study (see Fig. [3 compared to the predictions of the EM06 and 
EMIO models for analogous types (left and right panels, respectively). The prediction of the PLE model for total GNCs is shown in both 
panels for reference {grey dashed lines). Vertical dotted line: limiting magnitude for 100% completeness of our data. Region shaded in 
blue in the right panel: region where the predictions of the EMIO model are not valid for total /C-band GNCs. Consult the legend in 
each panel. [A colour version is available in the online version.] 
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Figure 10. Apparent colour-magnitude diagrams {B — K vs. K) for our observational sample in the three wide redshift bins considered 
in P13 (0.3 < z < 0.7, 0.7 < z < 1.1, and 1.1 < 2 : < 1.5). Red circles: red regular spheroid-dominated galaxies (E-SO’s). Green diamonds: 
red regular disc-dominated galaxies. Yellow stars: red highly distorted galaxies. Blue crosses: blue galaxies. Lines: theoretical trends 
followed by various galaxy types within the redshift interval of each frame, modelled assuming standard SFHs starting at Zf = 3 and 
characteristic physical pr operties according to observations. They have been obtained with the stellar population synthesis models by 
iBruzual Sz Charlotl ll2003l . see the text for details). Stellar masses of M* = 10^'^, 2 X 10^^, and 4 X IO^^Mq have been assumed for 
all models at 0.3 < z < 0.7, 0.7 < z < 1.1, and 1.1 < z < 1.5 respectively. Red thin solid: E-SO galaxy. Orange thick solid: Sa-Sab. 
Green dashed: Sb-Sbc. Blue dot-dashed: Sc-Scd. Purple three dots-dashed: Sd-Irr. Pink long dashed: Dust-reddened starburst galaxy 
with Ay = 3.0 mag. [A colour version is available in the online version.] 


red sequence, Sb-Irr’s are the blue galaxies in the data sam¬ 
ple, and the advanced stages of major mergers predicted by 
the model are compared to the galaxies with distorted mor¬ 
phologies identified on the red sequence. This shows that 
it is feasible to explain this assembly of massive E-SO’s at 


0.8 < 2 < 1.5 (revealed by GNC data) only through major 
mergers. 

We conclude that the slope change in total NIR GNCs 
is thus a relic of the assembly of a substantial fraction of 
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present-day massive E-SO’s at 0.8 < 2 < 1.5, in agreement 
with hierarchical theories of galaxy formation. 
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